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Abstract 



Quantum field tlieories and Matrix models have a far richer solution set 
than is normally considered, due to the many boundary conditions which must 
be set to specify a solution of the Schwinger-Dyson equations. The complete 
set of solutions of these equations is obtained by generalizing the path inte- 
gral to include sums over various inequivalent contours of integration in the 
complex plane. We discuss the importance of these exotic solutions. While 
naively the complex contours seem perverse, they are relevant to the study of 
theta vacua and critical phenomena. Furthermore, it can be shown that within 
certain phases of many theories, the physical vacuum does not correspond to 
an integration over a real contour. We discuss the solution set for the special 
case of one component zero dimensional scalar field theories, and make remarks 
about matrix models and higher dimensional field theories that will be devel- 
oped in more detail elsewhere. Even the zero dimensional examples have much 
structure, and show some analogues of phenomena which are usually attributed 
to the effects of taking a thermodynamic limit. 
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I. INTRODUCTION 



This paper is a consequence of our ongoing study of a new numerical approach to quan- 
tum field theory, known as the Source Galerkin technique. This method is based 
on systematically approximating the functional differential equations, or Schwinger- Dyson 
equations, satisfied by the generating functional in quantum field theory, while controlling 
the error with a weighted averaging procedure. While exploring this method, it came to our 
attention that it is necessary to carefully analyze the boundary conditions imposed on these 
equations. 

Here, we address this problem from a formal rather than a numerical point of view. 
For any theory with a large but finite number of degrees of freedom, there are a very large 
number of "theta" parameters characterizing solutions of the Euclidean Schwinger-Dyson 
equations, which are the quantum mechanical version of the equations of motion. Only 
one point in this parameter space corresponds to the usual path integral Q. It is therefore 
tempting to discard all the other solutions as unphysical. However, we shall argue that 
under many conditions the exotic solutions are physical. Symmetry breaking vacua provide 
the simplest examples of a situation in which exotic solutions are physical. Furthermore, 
critical phenomena have a natural interpretation in terms of the behavior of the full solution 
set in the thermodynamic limit. In this limit the solutions associated with many different 
boundary conditions coalesce. For other boundary conditions, the limit does not exist. 
Which boundary conditions have a thermodynamic limit and which ones coalesce depends 
on the theory's parameters in a manner which determines the phase diagram and the critical 
exponents. 

In section II we discuss the boundary conditions of some simple zero dimensional scalar 
theories. In addition to the usual integral solution of the Schwinger-Dyson equations, there 
are exotic solutions which are obtained from sums over various inequivalent complex contours 
of integration. Zero dimensional theories have many phases, all of which are continuously 
related. For symmetric actions these include symmetry breaking solutions. 

In section HI, we describe the boundary conditions of the Schwinger-Dyson equations 
for a general theory with a finite number of degrees of freedom. For a Euclidean field theory 
on a finite lattice without a boundary, the complete set of solutions may be obtained from 
different complex integral representations and have a simple correspondence to the classical 
solutions of the theory. The exotic solutions may be thought of as generalized lattice theta 
vacua. There is a still larger class of solutions which may be obtained on a lattice with 
a boundary, at which the Schwinger-Dyson equations may be modified. In the integral 
representation of the solution, the measure for fields at the boundary is essentially arbitrary, 
and may be set by surface terms in the action or by wave functions. However we will confine 



* Sometimes no point corresponds to the usual path integral. This is the case if the action is 
unbounded below. Such a situation occurs in Euclidean Einstein gravity and in certain matrix 
models 



1 



our discussion to generating functionals in the vacuum sector. These may be obtained on a 
periodic lattice, for which the exotic contours provide a complete set of boundary conditions. 
In the thermodynamic limit the real path integration with surface terms in the action gives 
some but not all of the solutions found using complex contours for a periodic lattice[|. There 
are compelling reasons to consider solutions represented by exotic contours, of which we list 
a few below. One reason is that certain aspects of critical phenomena are very naturally 
related to the behavior of the full set of solutions, including those which are unphysical, in 
the thermodynamic limit. The full solution set may only by obtained by allowing complex 
contours. Furthermore solutions corresponding to false vacua may be of physical interest 
but are not obtainable from the integration over real fields. In the case of matrix models it 
is essential to consider exotic contours to obtain even certain physical solutions with a real 
double scaling limit. This is true even when the matrix potential is bounded below for real 
eigenvalues. In this case there is no analogue of a spatial boundary at which the measure 
may be arbitrary. There are often many physical solutions having a real double scaling limit. 
One can obtain some of these solutions from the integration over real eigenvalues by making 
global perturbations of the action which are removed after taking the — > oo or double 
scaling limit. However there are also examples of solutions which can not be found this way 
01 m P • the other hand, the complete set of exotic contours does yield the full solution 
set There are situations in which the boundary conditions associated with a real 
solution in the double scaling limit are extremely unusual |^ 0. 

In section IV, we discuss phenomena which occur in the thermodynamic limit in the con- 
text of the Euclidean Schwinger-Dyson equations. In particular, we discuss the coalescence 
of solutions associated with different boundary conditions, and the absence of a thermody- 
namic limit for other boundary conditions. In some cases one can explicitly see how the 
solution set collapses, as in certain large matrix models. The solutions which survive in 
this limit do not necessarily coalesce with a solution associated with the integral over real 
fields, and sometimes the integral over real fields itself does not have a thermodynamic limit. 
In this paper we will not explicitly demonstrate the collapse of the solution set. Instead we 
will assume the collapse and show that it leads to several expected phenomena, such as 
the appearance of phase boundaries in the thermodynamic limit due to the accumulation of 
Lee- Yang zeroes. Critical exponents are determined by the manner in which the solution set 
collapses. We give a simple zero dimensional analogy in which tuning a coupling constant to 
zero shrinks the solution set and causes Lee- Yang zeroes to accumulate. Setting the coupling 
constant to zero in this example is analogous to taking the thermodynamic limit. We also 
show that the collapse of the solution set induces another set of well known equations which 
we will refer to as the Schwinger action principle!; 



^'Here we assume that for real fields the potential is not unbounded below. 

■^The Schwinger action principle can be used to produce all equations associated with the formula- 
tion of a quantum field theory including the Schwinger-Dyson equations. However, as a matter of 
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Here is a parameter of the the action 5", and Z is the partition function. If one only 
considers the usual solution of the Schwinger Dyson equation corresponding to the integral 
over real eigenvalues in a matrix model, or real fields in a field theory on a periodic lattice, 
then the action principle is obviously obeyed. However, as soon as one considers all of the 
solutions of the Schwinger- Dyson equations, it is no longer obvious that the action principle 
should be satisfied. For a finite system the action principle is independent of the Schwinger- 
Dyson equations. It is possible to satisfy the Schwinger-Dyson equations while violating the 
action principle. Since the solution set does not necessarily collapse to the usual integral 
solution, it is nontrivial to show that the action principle follows from the Schwinger-Dyson 
equations in the thermodynamic limit. 

In section V we discuss the effect of the collapse of the solution set on the form of 
the effective potential. For a finite number of degrees of freedom, there are an infinity of 
effective potentials, each of which has only single extremum. In the thermodynamic limit, 
it is possible to have a single effective potential with multiple extrema. 

II. BOUNDARY CONDITIONS OF SCHWINGER-DYSON EQUATIONS 

When there is no spatial boundary, the set of boundary conditions of the Schwinger- 
Dyson equations is determined by the behavior of the action at asymptotic values of the 
fields. The reasons for this will become clear shortly. For many theories a potential term 
dominates the action at large fields, and wc will begin by considering a theory with only a 
potential term. Let us find the space of solutions for a polynomial action with one degree 
of freedom, 

n -1 

S^Y.l9ict>' (2.1) 
1=1 

The generating function Z{J) of disconnected Green's functions satisfies the Schwinger- 
Dyson equation: 

12(91 gjiZT - J)Z{J) = (2.2) 

This is a linear differential equation of order n — 1, and has an n — 1 parameter set of 
solutions. One parameter corresponds to a trivial normalization factor, so there is an n — 2 



convenience, we will refer to equations exclusive of the Schwinger-Dyson equations as the Schwinger 
action principle or just the action principle. 
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parameter set of solutions with different Green's functions. There is a complete basis set of 
solutions with the integral representation 

Z{J) = d(t>e-^^'>'^+^^ (2.3) 

where F is a complex path such that 

/ d0_^e-^W+^^ = (2.4) 
Jr ad) 



d(f) 
or 

g-S(0)+J0 



The paths for which this is true is determined by the order of the polynomial potential. We 

require Re{gn(j)'^) — > +oo as 101 — > cxo, giving n domains in which the contour can run off 
to infinity (fig.l). One can then write down n — 1 independent contours which comprise a 
basis set. The contours may not be deformed into one another since e""^ is not analytic at 
\4>\ — > oo. For S = \g(i)^ + ... with g real, one choice of independent contours is 

F+ = [-00, 0] + [0, +ioo] 
F° = [-00, 0] + [0, +oo] 

r- = [-00, 0] + [0, -ioo] (2.6) 

with which we associate the generating functions Z"*", and Z~ . As a matter of definition, 
we will say that solutions at different values of g are in the same phase if they are related 
by the the Schwinger action principle, 

d Id* 

For the example above, a general solution of the Schwingcr-Dyson equations is of the form 
Z = aZ^ + bZ^ + cZ^ where a, b and c may be arbitrary functions of g. The solutions Z^, 
Z° and Z~ separately satisfy the action principle. Therefore the action principle is satisfied 
if a,b and c arc independent of g. A phase is labeled by the ratios of these three constants. 
However as Re{g) becomes negative, the contours F"*", F° and F~ are no longer convergent. 
For g real and positive, there are domains swept out by contours equivalent to F+, F° and 
F~. As one rotates the phase of g, these domains also rotate to maintain convergence. So to 
be more precise the action principle is satisfied if for infinitesimal variations of g there is a 
choice of contours which is held fixed. If one makes large changes in g, it may be necessary 
to change the contour. Thus the statement that two solutions are in the same phase if 
they arc related by the action principle makes sense only locally in the space of coupling 
constants. A large loop in coupling constant space can go around a phase boundary and onto 
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FIG. 1. The shaded area is the region of the complex <j) plane defined by cos(40) > 0, where 6 
is the argument of (j). Any path starting and ending at infinity within one of these four domains 
corresponds to a particular solution for the zero dimensional 0^ theory 



another Riemann sheet, thus changing the boundary condition. In the above example there 
is a phase boundary at 5^ = 0, and a 27i rotation of g forces a | rotation of the contour to 
maintain convergence. Note that in zero dimensions the Schwinger action principle is a very 
arbitrary way to define a phase. One could have required the generating function within a 
given phase to satisfy any other set of first order differential equations in the couplings. As 
we shall see later, the action principle is a consequence of the Schwinger-Dyson equations in 
the thermodynamic limit |. 

For the example of the action 5* = + the exact disconnected Green's functions 
(see Appendix I) satisfying both the Schwinger-Dyson equations and the action principle 
are given by; 

and 



sEven in zero dimensions, if one solves the Schwinger-Dyson equations in a weak coupling expan- 
sion, the action principle is automatically satisfied order by order. By performing a weak coupling 
expansion, one has chosen a restricted class of phases with a particular weak coupling asymptotics. 
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G2n+1 — A* 
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i2n!! 



i-t) 



n 



(2.9) 



vil,t) UiO,t) + pUiO,-t) 




where ^ = y f^' U{a.t) and V{a.t) arc parabolic cylinder functions defined in appendix 
I. The parameters a and p depend on the ratios of the coefficients, a, b and c, of the 
independent contours. By requiring Green's functions to be real, or that a be real and 
b — c*, one is still left with many exotic sohitions. These are labeled by the two real 
parameters a and p. Solutions with a ^ break the (f) —(f) symmetry of the action. Such 
solutions exist because the contours and break the symmetry. 

Since reality is not much of a constraint on the solution set, it might be tempting to 
require positivity of even Green's functions as well. This is a much stronger constraint, 
since only the solution corresponding to the real contour is manifestly positive. For a 
general boundary condition, one must do some work to ascertain positivity. However, it 
does not make sense to impose positivity before taking a thermodynamic limit, because the 
conditions for positivity may change in this limit. Our approach to the selection of boundary 
conditions will be to let the thermodynamic limit do most of the selection for us. We shall 
argue later that the solution set collapses in the thermodynamic limit. By this we mean 
that the solutions associated with certain continuous classes of boundary conditions coalesce, 
while other boundary conditions lead to no thermodynamic limit. The surviving solutions 
may then include some which violate reflection positivity (for Euclidean Schwinger-Dyson 
equations), and these may be discarded by hand. Note that sometimes complex solutions 
describe false vacua and are of physical interest. 

One might also try to constrain the solution set by requiring that the weak coupling limit 
smoothly approach the free field solution. This would mean that Z should have the formal 
representation 



This is a very strong constraint reflecting the fact that the Schwinger-Dyson equation be- 
comes flrst order at 5^ = 0. For the zero dimensional 0^ theory the equation is {pdj — J)Z = 
and, consequently, has only one solution up to a normalization. However, this constraint is 
ill deflned if there are any phase boundaries in g. In our zero dimensional example there 
is a branch point at = 0, so the result of taking a weak coupling limit depends on the 
path. One could take — > with the phase of g flxed. There is no physical reason why 
the partition function should have a free fleld limit as g ^ 0. In the symmetry breaking 
phase of a (70*^ theory, the Green's functions are singular at g = 0. For instance, the one 
point function at weak coupling is < >= \J~^- For our zero dimensional example, the 
requirement of a free fleld g {) limit with < would yield the symmetric boundary 
condition corresponding to an integral over the imaginary axis. 

It is enlightening to consider the weak coupling limit of the full set of solutions of the 
zero dimensional 0^ theory. There are some solutions to which the effective potential in the 
loop expansion is asymptotic. These correspond to sums over contours which pass through 
a single dominant saddle point. At each order in the loop expansion the effective potential 



Z{g,J) = e 




(0,J) 



(2.10) 
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has three extrema corresponding to each saddle point and a different class of boundary 
conditions 0. Among these are symmetry breaking solutions with the asymptotic expansion 
for the one point function given by 



<<P>=J-^{l + 0{g) + ...) (2.11) 




is one example of a solution with this expansion. Another example which is real is 
+ Z~ for g > and /i < 0. One can explicitly show that such a boundary condition 
leads to the symmetry breaking phase of a Tr0^ matrix model 0. It seems natural 

to conjecture that some generalization of such a boundary condition also corresponds to the 
symmetry breaking phase of a 0^ field theory. 

There is yet another class of solutions of the zero dimensional 0^ theory which does not 
correspond to the loop expansion. These arise when the contours cross degenerate dominant 
saddle points. These solutions are linear combinations of the ones associated with the loop 
expansion. Some of these are extremely singular at weak coupling. For instance for fi > 
and g > 0, the solution Z = Z^ + Z~ has a one point function with the expansion 



<0>=^^e+^(l + ...) (2.12) 

This is also the expansion for a similar solution in the case fi < 0, but the boundary 

condition is different, with contours rotated by |. The essential singularity e arises 
because < >= and the contribution of the two degenerate dominant saddle points 
at 6 = X — cancels in the denominator but not in the numerator. In one matrix models it 
is easy to see that boundary conditions of this singular type do not have a thermodynamic 
limit 0. However at this point of our discussion we already have some evidence that 
some the other exotic boundary certainly may be physical, the simplest example being that 
of the symmetry breaking solutions associated with the loop expansion. 

Our interest in 0^ theory is primarily as a toy model. The same analysis may be applied 
to obtain the solution set for non-polynomial actions, in which the order of the Schwinger- 
Dyson equations is not always obvious. A simple example is given by the loop equation 
of one plaquette QED, which is defined by the action S = —(3cos(j). The condition that 
vanishes as |0| — * oo restricts the domain of integration to satisfy Re{(3cos(f)) — * +oo 
asymptotically. For (3 real and positive, one choice of inequivalent contours which satisfy 
this requirement is given by 

Fi = [-icx), +ioo] (2.13) 



**The three extrema are a spurious feature of the loop expansion in zero dimension. In fact non- 
perturbatively there is a different effective potential for every boundary condition. This will be 
discussed in a later section. 
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and 



T2 = [-ioo, 0] + [0, 2n] + [27T, 27t + ioo] (2.14) 

to which we associate the partition functions Zi and Z2. Due to the periodicity in Re{(j)), it is 
easy to see that any other aUowed contour may be written as hnear combination of these two. 
The difference between Zi and Z2 is just the usual integral defined on the compact interval 
[0,27r]. Note that only this solution manifestly satisfies the condition that | < e**^ > | < 1. 
However along the lines of our previous arguments, we will not impose such conditions before 
taking a thermodynamic limit. 

Given that there are two independent contours, we expect that the Schwinger-Dyson 
equations should be second order. In this case, it is easy to construct the equations even 
though the action is not polynomial. By coupling sources J and J to the loop variables e*"^ 
and e"**^ one obtains the Schwinger-Dyson equations 

- dj) - 2{Jdj - Jdj)]Z{J, J) = (2.15) 

and the constraint 



djdjZ{J,J) = Z. (2.16) 
These equations may be rewritten in the form 

djZ = Z (2.17) 



P-2J 
and 

which, accounting for the normalization of Z, has a one dimensional space of solutions 
characterized by the above integral representations. Of course, in contrast to a theory with 
a local order parameter, one can not hope to learn very much about the phases of a gauge 
theory from a zero dimensional (one plaquette) model. Note that lattice gauge theories are 
examples of theories in which the large field behavior of action, and therefore the choice 
of boundary conditions, is not controlled by independently variable local terms, due to 
the constraints among the plaquettes. However the one plaquette model does illustrate 
something which is true generally. The "exotic" solutions of the loop equations for a gauge 
theory correspond to certain complexifications of the gauge group. Similarly, in hermitian 
matrix models the exotic solutions correspond to certain complexifications of the eigenvalues 
in the matrix integral, and in this case it can be proven that some of the exotic solutions 
are "physical" for certain values of the couplings, meaning that a real double scaling limit 
exists. 
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III. FROM ZERO TO HIGHER DIMENSION 



For theories in which independently variable local terms dominate the large field behavior 
of the action, it is straightforward to represent the higher dimensional theory in terms of 
a zero dimensional theory. Such representations are the starting point for both the strong 
coupling and mean field expansions. For example, consider the lattice Schwinger- Dyson 
equations of an interacting scalar field with action S = ^ Y.x,y Hx,y<t>x<t'y + Hx Vi.'t'x)'- 

(E i/x,,^ + V\^) - J^)Z[J] = 0. (3.1) 

A formal solution corresponding to the strong coupling expansion is given by 

Z[J] = ^-'y ^^^'^^ [] Zo( J^) (3.2) 

w 

where Zq[J) is any solution of the zero dimensional Schwinger- Dyson equation 

iy\^) - J)Z,{J) = 0. (3.3) 
The solution also may be written as 

Z[J] = [ l[dPxe-'^^^-""''^''UMJ- + P^)- (3.4) 

X X 

This may be thought of as a stochastically driven zero dimensional theory. A saddle point 
expansion in the auxiliary variable p yields a mean field expansion. We shall not make 
any use of these expansions here, but just mention them to illustrate the structure of the 
solution set in non-zero dimensions. In both the representations there is a product over zero 
dimensional solutions at each lattice site. It is not necessary to use the same integration 
contour at every lattice site, even if lattice symmetry is required of the solution. A sum 
over translations and rotations of solutions with inhomogeneous boundary conditions yields 
a solution which satisfies the lattice symmetries, and yet can not always be written in terms 
of an integral over a product of equivalent measures at each lattice site (fig.2). In a future 
publication we show that there are matrix model solutions with this type of boundary 
condition. In this case the lattice site is replaced with an eigenvalue label, and the integral 
representation does not involve a single measure for each eigenvalue 0. This means that 
the solution is not of the form, 



(3.5) 



where A„ are the eigenvalues, and A [A] is the Vandermonde determinant. If we consider 
only solutions of this form, which involve products of equivalent measures over the lattice 
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(a) 

(b) 

FIG. 2. A pictorial representation of translationally invariant solutions on a two site lattice, 
(a) is a product over equivalent measures at each lattice site, (b) may not be written as a product 
over equivalent measures at each site. 

sites, then the space of boundary conditions on a periodic lattice is the same as that of the 
zero dimensional theory. 

Note, however, that there is a class of theories for which there is no simple relation 
between the space of boundary conditions in a zero dimensional theory and a subspace of 
the boundary conditions on a periodic lattice. This class consists of theories in which the 
action at large values of the fields is not dominated by independently variable local terms. 
We have stated that gauge theories fall in this class. Another example of a theory in this 
class is the XY model, with lattice action 

S=-]^H,,yCOS{e,-ey). (3.6) 

For such theories it is usually much harder to de-numerate the boundary conditions of the 
Schwinger-Dyson equations, although there is a partial classification based on generalized 
lattice theta vacua which we discuss shortly. It seems likely that there is an equivalence or 
at least a large overlap between theories without a local order parameter, and theories for 
which the boundary conditions are not determined by independent local terms. Note that 
for a theory in which local terms control the boundary conditions, such as the 0^ theory, the 
behavior of the local order parameter < > is determined by the choice of these boundary 
conditions. If in the thermodynamic limit the solution set collapses, then certain boundary 
conditions and therefore certain values of the order parameter are selected. 

To argue that the exotic solutions involving complexified fields may in fact be physical, we 
have so far given the example of symmetry breaking. We will now try to make our arguments 
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more concrete. While the exotic solutions may appear very strange, many of them are 
actually familiar. For example, there is a simple relation between conventional theta vacua 
and exotic solutions with complexified fields, which we demonstrate below. For some simple 
cases it can be shown that a complete set of solutions of the Schwinger-Dyson equations 
may be generated from the complete set of classical solutions, both real and complex. This 
relation may be made explicit through the Borel re-summation of perturbative expansions 
about the classical solutions. For each classical solution, there is a loop expansion of the 
partition function of the form 



= T.i^J^re-'^''^n,h- (3-7) 

where i labels the classical solution, and is the number of degrees of freedom. This 
expansion is asymptotic, but its Borel transform 

B^ = J2 -^^nA"" (3-8) 

is believed to have a a non-zero radius of convergence under very general conditions in field 
theory. Term by term, the inverse of the Borel transform is given by 



^^e-* / e^B,{t) (3.9) 



hS'' 



However this is often ill defined, as Bi{t) may have singularities on the positive real axis 
due to instantons and renormalons. In some simple cases one can prove that as long as the 
contour of the t integration is closed, or begins at t = and ends at Re{t) = +oo avoiding 
all singularities, then one obtains an exact, and usually exotic solution of the Schwinger- 
Dyson equations satisfying the action principle^. Sometimes this can be simply proven by 
analytic continuation of the couplings from another situation in which the couplings are real 
and there are no singularities on the positive real axis. In Appendix A, we present a more 
direct and general proof for an arbitrary theory with one degree of freedom. Note that as 
long as the action has no flat directions, and the perturbative parameter is not infinite for 
any of the classical solutions, then the number of classical solutions is equal to the number 
of independent solutions of the Schwinger Dyson equations. As a simple example consider 
the theory S = + ^0^. Perturbation theory is an expansion in For non-zero /i, 
there are three classical solutions and three independent solutions of the Schwinger-Dyson 
equations. However at = there is only one classical solution = 0, but still three 
independent solutions of the Schwinger-Dyson equations. We shall address such issues as 



^^Note that sometimes the integral over t for an open contour is not always convergent [10|. We 
will not address this this subtlety here. 
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flat directions and renormalons in the context of boundary conditions for the Schwinger- 
Dyson equations in a later work. Assuming that perturbative parameter for any classical 
solution is finite, then an arbitrary solution of the Schwinger-Dyson equations consists of 
a linear combination of solutions generated from the classical solutions. Which solution is 
generated by a particular classical solution is, of course, somewhat arbitrary. For instance 
in the 0^ theory the resumed perturbation expansion about the classical solution = a/ ~ 



for /i > yields either or depending on how one avoids the singularity in B{i) due 
to the neighboring classical solution at t = S{(p = 0) = 0. 

The point we now wish to make is that an exotic solution of the Schwinger-Dyson equa- 
tions given by Z = CjZj, where the Zj are generated from classical solutions, is a generalized 
form of theta vacuum. The q play the role of theta parameters. If the theory has a space- 
time with a boundary, then these theta vacua are simply related to conventional theta vacua. 
In the usual approach a particular theta vacuum is selected by adjusting a surface term in the 
action and integrating over real fields. The surface term term effects the Schwinger-Dyson 
equations only at the space-time boundary, so in the infinite volume limit its role is only 
to set a boundary condition. It does this by putting a different weight on the contributions 
to Z coming from fiuctuations about different real classical solutions, which when Borel 
re-summed correspond to different complex contours. Thus, in the thermodynamic limit, it 
seems clear that the real integration with a surface term in the action is equivalent to an 
exotic solution Z — CiZi, although we shall not prove it with any rigor. 

There is a class of boundary conditions which we have not yet discussed. These boundary 
conditions arise because on a lattice with a boundary, the measure at the boundary is 
completely arbitrary. It docs not have to correspond to a surface term in the action which 
puts different weights on the expansion about different instantons, or equivalently to some 
choice of contours on a lattice without a boundary. For example in a field theory any matrix 
element of the form 



where the wave-functions ip and ip' are independent of J{x)^ satisfies the Schwinger-Dyson 
equations. These solutions arc obtainable from the real contour path integration with the 
wave functions setting the measure at the time boundaries. For instance in quantum me- 
chanics one would write 



Such boundary conditions are necessary to obtain generating functionals in the excited 
states, which can not be found from the generalized theta vacua discussed above. By con- 
sidering a free lattice field theory, it is easy to see that excited state solutions may not be 
obtained from exotic contours with periodic boundary conditions. Consider the Schwinger 
Dyson equations 





(3.10) 





(EC 



d d 



J.)Z[J] = 0. 



(3.12) 



y 



X 
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In Euclidean space on a periodic lattice (9^ is invertible and these equations have a unique 
solution up to a normalization. Unsurprisingly, there are no inequivalent contours in the 
integral representation. Therefore to get the solutions which in the continuum limit cor- 
respond to generating functionals in all the excited states, one must consider the same 
equations with a temporal boundary. However in Euclidean space, the excited state Green's 
functions are distinct from the vacuum Green's functions in that they do not satisfy cluster 
properties. For example < l|T(6(t)6(0))|l > behaves like | < l|O(0)|0 > pexp {Ei - Eo)t 
at large t. In this paper we will confine our discussion to generating functionals in the 
vacuum sector, for which the set of exotic contours on a periodic lattice is a complete set of 
boundary conditions. 

Note that different measures at a spatial boundary are not always available to set the 
boundary conditions. An example is given by the one matrix model. There are symmetry 
breaking solutions of the |TrM^ + . . . matrix model with real positive g which are not 
obtainable from the real contour by any deformation of the action [§[• However all 
the solutions may be obtained using the complete set of exotic contours 101 O- 

IV. SOLUTION SET IN THE THERMODYNAMIC LIMIT 

According to the previous arguments, any theory with a finite number of degrees of 
freedom will have multiple theta parameters which continuously relate all the solutions of 
the Schwinger-Dyson equations. In general the number of lattice theta parameters grows 
exponentially with the number of degrees of freedom. Thus it appears that there is an 
intractable number of solutions. Fortunately, not all of these solutions can survive, or 
are distinct, in the thermodynamic limit. This limit exists if the free energy is extensive, 
meaning that the partition function Z behaves like e'^^^ as the number of degrees of freedom 
or the volume of the system. A/", goes to infinity. In the thermodynamic limit, a linear 
combination of two solutions for which the free energy is extensive, Z = ae'-'^^^ + he~-'^^^, 
is either equivalent to the solution with the lowest free energy density, or does not have a 
thermodynamic limit if the real parts of /i and /2 are equal. Thus there are really only two 
solutions rather than a continuum of solutions labeled by the parameter |. In a Euclidean 
field theory the absence of a thermodynamic limit is reflected in the failure of the vacuum 
Greens functions to cluster. In an NxN matrix model it is reflected in the absence of a 
genus expansion in j^. 

The absence of this continuum of solutions is the reason that there is no spontaneous 
symmetry breaking in continuum quantum mechanics, even though on a finite discrete time 
lattice, the Schwinger-Dyson equations have real symmetry breaking solutions. Consider 
the partition function of the anharmonic oscillator, V{x) = |a;^ -|- with /x < 0, with 
periodic boundary conditions in time, and a period given by /3. There is symmetry breaking 
at every order in the perturbative expansion about x = \J^- Due to tunneling into the 
other well, the Borel transform of this expansion has a singularity on the positive real axis. 
By avoiding this singularity one obtains an exact but complex symmetry breaking solution 
of the Euclidean Schwinger Dyson equations. Furthermore since the free energy of this 
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solution is extensive in the perturbative expansion, its Borel re-summation is also extensive. 
In other words, for large /3 this solution is of the form Z = e~^^°, where Eo is complex. 
Since the Schwinger Dyson equations are linear, one can add other solutions in many ways to 
cancel the imaginary parts. For example Z + Z* is a real, symmetry breaking solution of the 
Schwinger Dyson equation. However the real parts of the free energy densities associated 
with Z and Z* are equal. Therefore the sum is not extensive, and the Euclidean greens 
functions do not satisfy cluster properties. In some cases it is possible to combine solutions 
with imaginary parts and still obtain an extensive solution. For example, it has been shown 



1^ that if one considers the instanton gas sum for the anharmonic oscillator, then the 
various complex parts arising from the attractive interactions of the instantons, and the 
Borel re-summations in different instanton sectors, all cancel, yielding a real solution which 
is extensive and symmetric. In this case a real and extensive solution is obtained by adding 
an infinite number of complex extensive solutions. There may be several ways to perform 
such infinite sums. In the quantum mechanical sine-gordon model, the different ways of 
doing so yield a periodic "theta" band, and in a field theory yield the conventional theta 



vacua. 

/>4 



Presumably in the broken phase of a field theory, there are real symmetry breaking 
sums over contours for which the thermodynamic limit exists, whereas for the integration 
over real fields there is no such limit. In this case one can still use the integration over 
real fields if one adds a symmetry breaking term at the space-time boundary, or a global 
symmetry breaking term which is removed after taking the thermodynamic limit. However, 
we wish to emphasize again that the full set of complex integration contours is sometimes 
necessary to account for all the solutions. This is especially clear in certain real double scaling 
solutions of matrix models @] [0 0, and also in false vacuum solutions of field theories. 
In the semi-classical limit, small perturbations of the action which are removed upon taking 
the thermodynamic limit tend to pick out quantum states which are are peaked about a 
particular choice among degenerate classical global minima. However, when there are non- 
degenerate minima, there are generally complex false vacuum solutions corresponding to the 
sub-dominant minima, and these may not be obtained from the real contour integration. 
As one evolves these solutions through a first order phase transition, they may be obtained 
from a sum over contours which is real. This discontinuous change in the list of boundary 
conditions which lead to a thermodynamic limit is typical of crossing a phase boundary. As 
we shall argue next, certain features of critical phenomena are also simply related to the 
existence of the generalized set of boundary conditions of the Schwinger- Dyson equations. 

In a one-matrix model is not difficult to see that in the thermodynamic limit, extensive 
solutions coalesce in such a way that there are either discrete sets of solutions or a countable 
number of theta parameters 0- It is difficult to prove this in a general lattice field 
theory. Instead we will assume that the solution set collapses in a way which depends on 
the coupling constants, and argue that this leads to expected physical consequences. One 
such consequence is the appearance of phase boundaries. In a finite system, as one makes 
large changes in the couplings one must change contours to maintain convergence. In the 
thermodynamic limit, one may have also have to make changes in the boundary conditions 
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so that this hmit continues to exist, and so that the solutions vary smoothly in the space of 
coupling constants. For a finite number of degrees of freedom this additional change in the 
boundary conditions (the sum over integral contours) violates the Schwinger action principle, 
as discussed in section II. However in the thermodynamic limit it is possible to make this 
change in the boundary conditions and satisfy the action principle if solutions with different 
boundary conditions coalesce. Let us assume that for every solution in the thermodynamic 
limit, there is a set of associated boundary conditions. Then a sufficient condition for two 
solutions at infinitesimally different values of a coupling to be related by the action principle 
is that these sets have some intersection. In order that a solution vary smoothly in the the 
thermodynamic limit, certain boundary conditions may become disallowed as one varies the 
couplings. However since there are generally large sets of boundary conditions associated 
with a solution, one can find boundary conditions which may be held fixed at least for small 
variations of the couplings. For large variations of the couplings one may be forced to change 
the boundary condition, but this does not violate the action principle the way it would for a 
finite number of degrees of freedom, in which case different boundary conditions correspond 
to distinct solutions. Thus by making a large loop in the complex plane of some coupling 
constant, one can return to a different boundary solution than the one with which one 
started^. Therefore because of the collapse of the solution set in the thermodynamic limit, 
phase boundaries can appear in the form of branch points in certain coupling constants. 
If one were to keep the boundary conditions fixed globally, then there would have to be a 
barrier to smooth analytic continuation, which according to the standard lore is formed by 
the accumulation of Lee- Yang zeroes P] of the partition function. 

To clarify the arguments above, we describe a simple zero dimensional analogue, in 
which the solution set shrinks accompanied by the accumulation of Lee- Yang zeroes and the 
appearance of a phase boundary. Consider a solution of the Schwinger-Dyson equation for 
the action S{(j)) = -|- + . . .. If the action principle is satisfied, then the generating 
function Z is analytic in A (except at infinity) as long as g does not vanish. There is no need to 
rotate contours to maintain convergence as one changes the phase of A. However if one sets g 
to zero, then the order of the Schwinger-Dyson equation {gdj + Xdj+^dJ—J)Z = drops by 
one, and the space of solutions becomes smaller. In this sense, the limit g ^ is analogous to 
the thermodynamic limit. Different boundary conditions give different asymptotic behavior 
for small g. Some solutions do not exist in the g 0'^ limit, while others coalesce. The 
sets of boundary conditions which survive and coalesce in this limit depends on the value 
of A. Consider the limit g ^ with g real and positive. The sets of equivalent convergent 
integration contours for g ^ differ from those for (7 = (fig.3). The intersection of these 
sets determines how solutions of the theory behave as g is set to zero. If the set of equivalent 
contours associated with a solution for 7^ has no intersection with the set of equivalent 
contours associated with any solution for g = 0, then there is no — limit. For instance 



^-l- Starting from a physical solution, such analytic continuations may lead to solutions which do 
not satisfy physical constraints, or to other physical solutions. 
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the integral over the real axis behaves as Z ~ exp{^^) as g ^ 0^ with A real and positive. 
For certain values of A there are solutions for which the intersection consists of contours 
which can be closed in the A0^ theory, giving Z — * as (7 — > 0. Therefore at these values of 
A, there are many solutions which coalesce as g ^ 0^ due to the linearity of the Schwinger- 
Dyson equation. At g = maintaining convergence forces one to rotate the contours as 
one rotates the phase of A by 2tt, and the solutions become triply sheeted in A with a 
phase boundary at A = 0. As one changes the phase of A at small but non-zero g, small 
violations of the action principle allow one to move among solutions in a way which at g = 
corresponds to analytic continuation among Riemann sheets in the complex A plane. This 
amounts to altering the boundary conditions to prevent Stokes phenomena from occuring 
on the positive real g axis as one changes the phase of A. If one were to keep the boundary 
conditions (contours) fixed at non-zero g, then Z would be a single valued function of A and 
there would have to be some obstruction to analytic continuation of A in a small g expansion. 
One way this obstruction can appear is by the accumulation of Lee- Yang zeroes of Z along 
a Stokes line in the complex A plane as g 0^. At small but non-zero g one can change the 
location of the Stokes line by changing A or by changing the boundary conditions. By the 
appropriate change in the boundary conditions as one varies A, one can prevent the Stokes 
line from crossing the positive real g axis. The reason it is possible to make such changes in 
the boundary conditions with only small violations of the action principle is that solutions 
with different boundary conditions coalesce as g ^ 0. 

For the example we have given, it is difficult to see if the accumulation of zeroes explicitly. 
However there is an even simpler example. Consider starting with the theory S = -^^^ + 
and take the limit A — > 0"*". In this case the partition function is exactly given by 

Z(A,/i) = X-h~^{aM-^) + (3Bi{-^)) (4.1) 

4A3 4A3 

where Ai{x) and Bi are the two independent Airy functions. The combinations which 
survive in the A ^ 0"'' limit approach the free field solution 

Z(0,/i) = ^ (4.2) 

At finite A there is no branch point in /i, although the ^JJi analytic structure does appear 
at any finite order in a large A expansion. One way this can occur is if Z has an infinite 
number of zeroes in A which become dense at infinity. As A — > 0"^ we expect these zeroes 
to accumulate also at finite values of A. Because of the in the argument of the airy 
function, it is clear that this is indeed the case, with zeroes accumulating on the negative 
real axi^ As one rotates the phase of /i, the lines of zeroes rotate in the complex A plane, 
unless one simultaneously changes the boundary conditions, namely a and (3. 



^Of course this is a somewhat trivial example of a phase boundary, since the difference between 
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(a) (b) 

FIG. 3. (a) The allowed angular domains of integration as ^ oo are drawn for S = gcj)'^ + . . . 
in white and for S = Xcfy^ + ■ ■ ■ in black, where g and A are real and positive. The intersection of 
these domains determines the behavior of solutions as g ^ {). (b) is the same diagram drawn for 
a complex A. The solution associated with the depicted contour has a vanishing partition function 
as 5 — > 0, keeping the phase of g fixed. 

It is quite generally true that the collapse of the solution set of the Schwinger- Dyson 
equations is intimately connected to the accumulation of Lee- Yang zeroes and the appearance 
of phase boundaries, whether this collapse is brought about by tuning parameters to special 
values, or by taking the thermodynamic limit. The collapse of the solution set forces one 
to change the boundary conditions while making large variations of the couplings in order 
that solutions depend smoothly on the coupling constants. At the same time the collapse 
allows one to make these changes without violating the action principle. Consequently it 
becomes possible to have solutions with branch point singularities, and critical exponents 
are determined by the manner in which the solution set collapses. It can also be shown in 
certain matrix models that analytic continuation of a solution with a real double scaling limit 
and very conventional boundary conditions (real axis integration) leads to other solutions 
which also have a real double scaling limit, but very unconventional boundary conditions 
0. These unconventional boundary conditions can not be written by assigning a single 
equivalent sum over integral contours to every eigenvalue. In other words, these solutions 
may not be written in the manner of figure 2(a). 

We have argued above that the collapse of the solution set in various limits forces one 
to make changes in the boundary conditions when one varies the couplings, and that these 
changes may be consistent with the action principle. While this is clearly true in the zero 
dimensional examples, the arguments that it is true in a thermodynamic limit were only 



the two Riemann sheets in /x is only a normalization. A 2tt change in the phase of /i leads to a 
rotation of the contour by vr, which changes only the sign of Z and not any of the Green's functions, 
which have a pole rather than a branch point at /x = 0. 
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heuristic. Below we show that the action principle follows from the Schwinger Dyson equa- 
tions in the thermodynamic limit. For a finite number of degrees of freedom the action 
principle and the Schwinger Dyson equations are independent. The manner in which a 
solution changes as one varies a coupling is arbitrary for a finite number of degrees of free- 
dom, since at any fixed coupling all the solutions are continuously related. However let 
us assume that in the thermodynamic limit there are only discrete sets of solutions of the 
Schwinger-Dyson equations at fixed values of the couplings, and none that are continuously 
related. Let us also assume that solutions vary smoothly as one varies the couplings, with 
exceptional points corresponding to phase boundaries. We shall define the couplings ga so 
that they appear holomorphically in the action and the Schwinger-Dyson equations. The 



solutions will then be holomorphic functions of the couplingsQj. Therefore some first order 
differential equations of the form 

d^Z = {^-k^)Z[g,J] = 0. (4.3) 

are automatically induced when the solution set becomes discrete^. More precisely we 
should say that there is a discrete set of solutions for / = jjlnZ as A/" ^ cxd, and it is 
this quantity which satisfies some set of first order differential equations in the couplings. 
The solutions for Z associated with one of the discrete solutions for / may be written as 
Z = e~^^ + b where 5 is sub dominant in the large M limit. There are generally many 
possibilities for b. Therefore strictly speaking the Z's associated with a particular / in 
the M ^ oo limit will satisfy an equation of the form ([4.3|) up to some arbitrary terms 
subdominant to e~^^ . These arbitrary terms depend on how one chooses 5, given /. 

It remains to show that these equations correspond to the action principle. We shall 
assume is a linear operator, since a generic nonlinear term would be either dominant 
or sub-dominant to ^-Z^ in the thermodynamic limit. For instance Z"^ ~ e~'^^ while 

■Tp-Z ~ e~-^^ , as TV — >■ cxD. Thus we have 

dga ' 

Let us denote the Schwinger-Dyson equations by LiZ = 0, where the index i labels the 
degrees of freedom. Since Oa is linear, [Oa,Li]Z = 0, up to terms subdominant to e~-'^-^. 
We decompose Oa into a piece which commutes with all the Lj and a piece which does not; 



*** While analyticity holds in the thermodynamic limit, it does not imply that analytic continuation 
of a physical solution around a phase boundary yields a solution which is still physical 0|. 

tt^If the solutions in the thermodynamic limit are not discrete, but have a countable number of 
continuous theta parameters, then the Schwinger action principle is not induced, but may certainly 
be consistently imposed. 



18 



(4.5) 



where 

[Ha,Li]^0 (4.6) 
and the derivative with respect to the couphng is included in Ha] 



Ha=^-... (4.7) 
dga 



With this decomposition one has 



lim e+^f[da, U]Z = lim e+^^L.MaZ = (4.8) 

AT— »oo N—*oo 

Therefore, Af„Z may be written as the sum of any solution of the Schwinger Dyson equation 
Z' and a term subdominant to e"-'^-^. However for the equation OaZ = to make sense in 
the thermodynamic limit, Z' can not be dominant with respect to Z in the J\f ^ oo limit. 
Therefore Z must be an an approximate eigenfunction of M^, 

MaZ - ^{gk)Z (4.9) 

up to terms subdominant to e~^^ , which means that 

{Ha + l{gk))Z = Q, (4.10) 

up to terms subdominant to e'-'^^ . 

The set of operators Ha of the form -tt- ... which commute with Lj are in fact the 

operators associated with the action principle. To see this, consider an arbitrary zero dimen- 
sional theory with the action S = J2a9afa{4>)- The operator associated with the Schwinger 
Dyson equations is 

L = E9aadj)-J (4.11) 

a 

and the operators associated with the action principle are 

W,+M9.) (4.12) 

These operators commute. For instance with S — one finds that 

[§-g + \dj,gd'j-J]^0 (4.13) 

This commutation relation generalizes in an obvious way to theories with more than one 
degree of freedom. Note that one can add a term which is constant in the source J, but 
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dependent on the couplings, to any of the operators associated with the action principle 
without changing the commutation relations with L. There are some restrictions on what 
this additional term may be, since the change Ha —>■ Ha + 7a(5'i, 5*2, •••) leads to 

[H^,H^] = -^^--^^ (4.14) 

which must vanish if HaZ = The addition of a term 7^ satisfying this constraint has no 
effect on the action principle, since 7 may be absorbed by a change in the normalization of 
Z 

Z ^ e-/''^^"^"(^')Z (4.15) 

If the Schwinger action principle is written in terms of normalized Green's functions, the 7 
term does not appear at all. Thus the action principle follows from the Schwinger Dyson 
equations in the thermodynamic limit. 



V. ANALYTIC STRUCTURE OF THE EFFECTIVE POTENTIAL 

In this section we discuss the effect of the collapse of the solution set on the effective 
potential. In general, for a finite number of degrees of freedom, there is a continuous 
set of effective potentials, with extrema associated with a single or a small discrete set 
of solutions of the Schwinger-Dyson equations. Recall that the effective potential may be 
defined as follows. Take for example the (f)^ theory with action S — ^(p^ + ^cp"^ . In terms of 
0(J) = ^lnZ(J) the Schwinger-Dyson equations may be written as 

J = X(P\J) + li(P{J) + Xdj(P{J) (5.1) 

which has a one parameter class of solutions. Given a solution 0( J) the associated effective 
potential r(0) is then defined by the solving the equation 

j = r'(0). (5.2) 

This is an algebraic equation rather than a differential one, and in general r(0) can not yield 
the full continuous one parameter class of solutions, but only a discrete subset. There is an 
extremum of r(0) associated with each branch of 0(</) in the complex J plane. However in 
zero dimensions Z is analytic in J, and (j){J) will in general have an infinite tower of poles 
in J. In our example these poles arise from the infinite tower of zeroes of 

Z{J) = X-h-^^^^'M{^ + ^) (5.3) 

4A3 A3 

in the complex J plane. Since there are an infinite number of poles but no branch points 
in J, the analytic structure of r(0) is somewhat odd. Because of the infinite number of 



20 



poles, = cxD corresponds to an infinite number of discrete values of J = T'{(j)), so T{(j)) 
has an infinite number of Riemann sheets in the complex plane. Also, since 0(</) has no 
branch points in J, there is a unique value of (f) at which J = r'(0) = 0. Therefore there 
is a different effective potential for every solution. In general one must Legendre transform 
with respect to a variable in which the partition function is multiply sheeted in order to 
get an effective potential with multiple extrema. Note that at any finite order in the loop 
expansion, one gets an effective potential which appears to have more than one extremum. 
In zero dimensions this is a spurious feature of the loop expansion, which may be recast as 
an asymptotic expansion for large J. For instance, in the A^^ example, the loop expansion 
corresponds to writing 

J = X(f)\j) + J) + nXdj(f){J) (5.4) 

Perturbing in powers of h, and setting h = 1 gives an expansion of the form 

(p{J) = a{J - Jc)^ +b + c{J - Jc)^ + ... (5.5) 

where a,b,c. . . and Jc are functions of /i and A. At any finite order there is a square root 
branch point at J = Jc. In a thermodynamic limit however, it is possible that the analytic 
structure apparent in the loop expansion also corresponds to the analytic structure of the 
exact solution. In a limit in which the Lee- Yang zeroes of Z( J) accumulate to give a multiply 
sheeted solution, a single effective potential will describe several solutions. Furthermore, if 
solutions coalesce as they are expected too, then so will the associated effective potentials. 
There will not be a continuous class of effective potentials, except in cases in which a theta 
parameter survives. 

It is generally difficult to find the boundary conditions associated with a certain solution. 
However for the sake of illustration, let us consider the three extrema of the effective potential 
in the broken phase of 0"^ field theory, and make an educated guess as to the boundary 
conditions associated with each of these extrema. The two symmetry breaking extrema are 
probably obtained in the continuum limit of the lattice theory using contours such as 

E[_d<l>x + Yl[ (5-6) 



n / d<p, (5-7) 
V""'r-+r+ 



or 



where F"^ and F~ are as defined previously in ( p.6|) . 

The effective potential has another extremum at = 0, corresponding to an unphysical 
solution which survives in the thermodynamic limit. It is somewhat more difficult to guess 
which boundary conditions this solution corresponds to. We need another boundary con- 
dition which has a thermodynamic limit for n < and which is symmetric under — 
The lattice theory has a simple property which helps to find this boundary condition. The 
lattice action is 
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^ = 2 E Ky+i^<PAy + E + - DW (5.8) 

x,y,[i. X 

Where D is the dimension. Consider an arbitrary boundary condition and rotate the con- 
tours by ^ and ^ at alternate lattice sites. This is equivalent to keeping the contours 
fixed, but flipping the sign of the cjp' term and rotating the phase of the source term by 
i and —i at alternate lattice sites. The interaction and hopping terms are invariant. When 
the sign of the 0^ term is flipped fi — > —[i + 2D, so this mapping relates the broken and 
unbroken phases. Note that this is not a duality in the usual sense, since the boundary con- 
ditions are changed as well as the parameters of the action. If for < we choose contours 
of integration which get mapped to the real contour, then via this mapping the solution 
is equivalent to the one obtained from the real contour in an unbroken phase. The only 
difference is in factors of —1 for certain Green's functions, due to the rotation of the phase 
of the source term. We conjecture this solution corresponds to the local maximum of the 
effective potential. It would certainly not be surprising, since for // < in zero dimensions 
the local maximum of the effective potential at any finite order in the in the loop expansion 
corresponds to the solution obtained by integrating along the imaginary axis. 



VI. CONCLUSIONS 



The phase structure of bosonic field theories and matrix models appears to have a very 
natural interpretation in terms of the multiple boundary conditions of Schwinger- Dyson 
equations. In the thermodynamic limit, phenomena such as the accumulation of Lee- Yang 
zeroes and the collapse of the solution set are intimately related. There are several interesting 
questions which we are pursuing, of which we will list only a few. It is not known whether or 
how these ideas may be generalized to fermionic theories. The Schwinger- Dyson equations 
of a fermionic theory may be written, as in the bosonic case, as a set of recursion relations 
among Green's functions. However for fermions on a finite lattice there are countable number 
of Green's functions. Therefore these recursion relations must truncate. Because of this 
truncation, there is not a large set of solutions when the number of degrees of freedom is 
finite. This may have interesting consequences for the large field behavior of the action in any 
bosonized version of the lattice theory. Another question concerns the relation between the 
various boundary conditions and the phases of a theory with a nonlocal order parameter. 
As yet we have not been able to construct such a relation, though we conjecture that it 
exists. We are also pursuing the question of precisely how renormalons are related to the 
multiple boundary conditions of the Schwinger- Dyson equations. In the simple examples 
considered in this paper the perturbative parameter is finite, and there is a one to one 
relation between instantons and a basis set of solutions of the Schwinger-Dyson equations, 
with the only singularities in the Borel transform arising from instantons. When there are 
renormalon singularities as well, then there are more inequivalent ways to avoid singularities 
when integrating over the Borel variable. This leads to the conjecture that there are more 
boundary conditions for Schwinger-Dyson equations than instanton counting would indicate. 
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VII. APPENDIX I 

A. Exact solutions of zero dimensional (p^ theory 

We have discovered, not surprisingly, that considerable effort has been devoted to solving 
the zero dimensional (p'^ theory) [0, [|T6[, [|T^, |18|, ||T9|, |20|. Much of what we present in 



this section has been known in one form or another. It is the intent of our presentation to 
bring new clarity and completeness to the solutions of this model in a form consistent with 
the point of view of this paper. Here we construct the exact Green's functions for the full 
solution set of the zero dimensional theory with action S = + ^0^. The normalized 
disconnected Green's functions G„ are defined by the Taylor expansion of the generating 
function, 

^ = ^G^J"" (7.1) 
The Schwinger-Dyson equations may be written as the recursion relation 

gGn+3 + fJ'Gn+i - nGn^i = (7.2) 

The initial conditions for this recursion relation are the one and two point function. The 
Schwinger action principle may be written as 

-;GiGn + -G'n+4 = (7.3) 

og 4 4 

and 

^ - Ig^G^ + = (7.4) 

on 2 2 

Since the Schwinger-Dyson equations allow one to write all Green's functions in terms of 
G*i and G21 the action principle may be rewritten as a closed set of first order differential 
equations for Gi and G2-, which turn out to be soluble in terms of parabolic cylinder functions. 
Instead of working with Gi and G2, it is convenient to define the ratios 
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for n even, and 



<5n = (7.6) 

for n odd. Wc shall set = 1 for convenience. It can be reintroduced later via the Schwinger 
action principle for or equivalently by taking ^ to be dimensionless and (j) to have 

dimensions of In terms of the it!„, the Schwinger-Dyson equations read 

(2n + l)gRnRn-i + Rn-i -1 = (7.7) 

A combination of the Schwinger action principle and the Schwinger Dyson equations yields 
the Ricatti equation 

+ (1 + 2g)K + (2n + l)gRl -1 = (7.8) 



which by the change of variables, 

2t d 



Rr, 



-^/n(ei*V(t)) (7.9) 



2n + 

where t = leads to the defining equation of the parabolic cylinder functions; 

nt)-{\t' + n)m^O (7.10) 

This equation has two independent real solutions f/(n, t) and V{n.t). We give the definitions 
and some useful properties of these functions in the next section. One may also use U (n, t) 
and U{n, —t) as an independent set. Therefore we have 

fit) = anU{n, t) + PnU{n, -t) (7.11) 

which yields 

-anU{n+l,t)+(3nU{n+l,-t) 

ar,U{n,t)+f3r,U{n,-t) ^ ' ^ 

The Schwinger-Dyson equations relating Rn to Rn+i constrain the coefficients and 
so that we get 

o _ U{n+l,t) + {-ir+'pU{n + l,-t) 

U{n,t) + {-l)^pU{n,-t) ^ ■ ^ 
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The parameter p is not fixed and sets one of the two boundary conditions of the Schwinger- 

Dyson equations. 

Since Go = 1, the even Green's functions are given by 

G,^ - (2^ - 1)!!^ um+pU{n,-t) ^^'^^^ 
One can compute Qn in the same way we have computed Rn giving, 

- V(n - I, t) + i-l)nu;^V{n - 1, t) ^'-^'^ 

However in this case uj is not a free parameter. The Schwinger-Dyson equation gG^+Gi = 0, 
gives Qi = — ^ = — t^, and a; = cxd, so that 

nV[n - ^,t) 

To get the odd Green's functions from Qn one also needs to know Gi. Due to the Schwinger 
action principle, 

f (7.17) 

Using G5 = 8Q1Q2G1 and the above solutions for the even Green's functions, we find that 
<'«' = f(,+ l) + l£^iL4z4ik4)G. (7.18) 



Writing 



One gets. 



so that 



dt V 2U{0,t)+pU{0,-t)^ 



u{o,t) + pUio,-t) ^^-^^^ 



Thus the odd Green's functions are given by 



_2n!! ..n Vin + jt) t^eV 
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Since these are normalized Green's functions, the values of a and p are related to the ratios 
of the coefficients of the integral contours. It is amusing to note that there are classes 
of contours among which the even Green's functions are the same, and the odd Green's 
functions differ by an overall normalization. These are obtained by varying a with p fixed. 
This is a refiection of the fact that for the symmetric action, the Schwinger Dyson equations 
do not couple the even with the odd Green's functions. The coupling comes only from the 
action principle for the odd Green's functions, 

d 1 1 

T^G'2n+l ~ —G^G2n+\ + -G2n^^ = (7.23) 

og 4 4 

One can change the odd Green's functions by a global factor without effecting the above 
equations or the linear Schwinger- Dyson equations for the odd Green's functions. 



B. Truncated series and the path integral solution 



It has already been noted by several authors ( |T^) that solving the zero dimensional 
0^ Schwinger-Dyson equations 

by a truncated Taylor series leads to a solution Z'^ with smooth limit (7 — > 0"*" for > 0. 
It is straightforward to check that truncation causes odd coefficients in the series vanish. 
From (|7.7]), one obtains a continued fraction expansion for the ratios Rk 

Rk = . p A; = l,2...,iV-l , R^ = , N ^ 00 (7.25) 

1 + (2/c + l)gRk+i 



Rn = 7K (7.26) 

. , i2n + 3)g 
^ {2n + 5)g 

A truncated series provides then a rational (Fade) approximation for the Green's functions, 
determined by the classical convergents of the continued fraction ( |7.26| ). By Fincherle's 
theorem, (|7.26|) converges to the minimal solution of ( [7.2|) , the unique solution G^ 7^ such 
that ^ 

m—>oo (j' 



for any other linearly independent solution G'^. In fact, one can also prove |]14[ that ( |7.26| ) 
converges to a Stieltjes function: a function F{z) admitting the following integral represen- 
tation 
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F{z) = I dt-^^c^f^a^i-zr , ar,= rdtty{t) (7.27) 
Jo 1 + zt ^ Jo 

with pit) > Vt G [0, oo] . The series is asymptotic for z ^ and — vr < arg(z) < vr and 
the coefficients cr„ are moments of a positive density function p[t). Moreover, F[z) is the 
Borel sum of the asymptotic series if this sum exists. 
It is straightforward to prove that 

that is, equation ( [7.13| ) with p = 0. Note that Rn G [0, 1], where the maximum occurs at 
(7 = 0. Taking into account the relationship of the Green's functions with the ratios Rn-, and 
using the integral representations of the parabolic functions given in the next section, one 
finds that the truncated solution Z'^ converges to the usual path integral solution 

= — — --- (7-29) 



> e 



which, as stated, is the unique solution of ( [7.24|) with smooth limit g ^ 0^ for fi > 



C. Parabolic Cylinder Functions 

The parabolic cylinder functions U{a,t) and V{a,t) are independent (real) solutions to 
the equation 

^-(lt2 + a)y = (7.30) 

These function satisfy 



la 1,2 ,1 lllo. 



U{a, t) = 2-3- 2 e-4* ^(-a + -, -, -t') (7.31) 



7iV{a,t) =T{a + ^){simiaU{a,t) + U{a,-t)) (7.32) 



where is a hyper-geometric function defined in [^, pp. 503-535. 

1 1 19 1 1 

U(a, t) = — / e-*"-2« s^-^ ds a> — (7.33) 

r(a + ^) Jo 2 

The change of variables s = cfP' /2t transforms ( |7.33| ) into 
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Properties of these functions are listed, for example, in [|13[, pp. 685-720, and [|1J]. Some 
useful relations are 

U\a, t) + U{a, t) + (a + ^)U{a + 1, t) = (7.35) 

V'{a, t) + V{a, t) - V{a + 1, t) = (7.36) 

f/(a + 2,t) = ^^(^il^f^(^±M (7.37) 



22"+4r(ia + I) 



When is a non negative integer 



2 Ci j-2 



V{k + ^,t) = Hlit) = e-* ^e* 

(7.39) 

1 T 2 poo 

U{k + -,t) = -e- dz{t - zf e-^'^'' (7.40) 



2' ' k\ 
Asymptotic expansions for large t 

for |arg(t)| < ^ir 

Via,t) ^ Al^S^~Hl + + + . . .) (7.42) 

^ ' ' \7i ^ 2t^ 24^4 ' ^ ' 

for |arg(t)| < |7r 

VIII. APPENDIX II 

A. Borel re-summation and exotic solutions 

Here we demonstrate the relation between various Borel re-summations and the exotic 
solutions of the Schwinger- Dyson equations for an arbitrary polynomial action in a zero 
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dimensions, S{(j)) = ^0"'. It is shown that the by a suitable choice of integration contour 
in the Borel variable, one obtains an exact solution of the Schwinger-Dyson equation which 
satisfies the Schwinger action principle. The loop expansion about the the classical solution 
0Q, yields the following contribution to the generating function; 



We shall assume that all the coefficients c^are finite. This precludes actions with flat direc- 
tions, 5"' 7^ 0. Note that one obtains the same series starting from any integration contour 
for which has constant phase, and which passes through (p = (p^ but no other classical 
solution of equal or lower action. This series is asymptotic, but its Borel transform defined 
by 



has a finite radius of convergence. In the zero dimensional theory it converges to 

where in the vicinity of t = the contour C encloses, in the opposite sense, the two poles 
4>a,i{gk,t) which coalesce to (paidk) at t = 0. All the other poles are taken to lie outside 
the contour. The Borel transform has a singularity when one of the exterior poles coalesces 
with one of the interior poles, which occurs when t is equal to the action of a neighboring 
classical solution, t = S{4>a')- Doing the integral, we may also write. 

Bait) = E (-1)^^77T^ (8-4) 

where (f)a,i{t) are the two poles which coalesce at t = ; 0o,i(O) = (pa- 



Thus far everything we have said is standard We now exhibit an exact relation 

between the Borel re-summations and the exotic solutions of the Schwinger-Dyson equations. 
We invert the Borel re-summation by writing 

Z. = e-^^(^'^) Jdte~r.Ml (8.5) 

with an as yet unspecified integration contour in the complex t plane. We set h = 1 in what 
follows. If Za satisfies both the Schwinger-Dyson equations and the action principle then it 
is annihilated by the operators; 

L = J2(^- ^)9ndg„_, - gi (8.6) 

n=2 
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and 

dgn n dgl 



9 19" , , 

Hn^^--^ (8.7) 



It is convenient to define the quantity 

F« = / dte-'^ (8.8) 

If LZa = then CF^ = where 

i:=^(n-l)^„D,„_,-^i (8.9) 

n=2 

and 

Before proceeding, we list several simple but useful identities. Due to the equation of motion, 
S'{4>a{gk)) = 0, one has 

^ S{^a) = -C (8-11) 



OQn n 
Another obvious identity is 

Ki{t^^^9k)^4>a{gk) (8.12) 

Two other identities are found by differentiating the relation t — S{(f)i{t, gk)) — S{4>{gk)) with 
respect to t and g^- 



^-<.(t..>)= °'^£^-!r^?^^-" (8.14) 



The equations of motion are used again in deriving the last equation. Using these identities 
the quantity 



CF^I dte-\-iy-[Y,{n - l)gnDg„_, - gi]cf>a,i{t) (8-15) 



d 

n=2 

may be rewritten, after much algebra, as 
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/ 



-t 



(-l)T5nC-^|0a.(t)] 



(8.16) 



which vanishes for any contour in t via the equations of motion I^nS'n^a ~ 0. Note 
that it was also not necessary to encircle the zeroes of i — <S'(0) + S{4>a) in any particular 
way when doing the initial integral. The factor of —1* was certainly not required to 
satisfy the equation LZ = 0. Note however that this equation is some combination of the 
Schwinger-Dyson equation and the action principle. We must also see if the action principle 
is separately satisfied. To this end, consider the quantity 

[klD,p,, + {k + 1)D,,JF = / dte-'- Y.{-mklD,,D,, + {k + l)D,,Jcl>^,,{t) (8.17) 



which will vanish if the action principle is satisfied. Using the same identities this quantity 
may rewritten as 



/ 



dt 



d_ 
di 



_d__d_ 
dgk dgi 



E(-i) 



'e-^kl 



d d 
dgk dgi 



.18) 



<9r 



which vanishes for several choices of contours in the t plane. The contour may be closed 
around some Borel singularities, or it may begin and end at Re{t) = +oo winding around 
some singularities at finite t. Alternatively the contour may begin ai t — and end at 
Re{t) = +00. The latter choice is viable because at t = the two poles 4>a,i{t) coalesce 
and the factor of I]j(— 1)* then causes the t = boundary term to vanish. As long as some 
such contour is used, then one obtains an exact solution of the Schwinger-Dyson equation 
satisfying the action principle. We suspect that the analysis above generalizes to theories 
with a greater number degrees of freedom, provided there are no flat directions, though as 
yet we have not constructed a proof. 
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